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Analytical expressions for the entanglement measures concurrence, i-concurrence and 3-tangle in 
terms of spin correlation functions are derived using general symmetries of the quantum spin system. 
These relations are exploited for the one-dimensional XXZ-model, in particular the concurrence and 
the critical temperature for disentanglement are calculated for finite systems with up to six qubits. A 
recent NMR quantum error correction experiment is analyzed within the framework of the proposed 
theoretical approach. 
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I. INTRODUCTION 

Quantum entanglement was already pointed out by 
Schrodinger pj to be a crucial element of quantum me- 
chanics. Research was refocused on quantum entangle- 
ment in the last fifteen years because the field of quantum 
information theory (cf. @, Q) developed rather quickly. 
Recent papers concerning entanglement in quantum spin 
systems address questions about the maximum entan- 
glement of nearest neighbor qubits belonging to a ring 
of N qubits in a translationally invariant quantum state 
Q, the dependence of entanglement between two spins 
on temperature, external magnetic field strength and/or 
anisotropy for the one-dimensional isotropic Heisenberg 
model ISlEllJlHinj, Isin § model EL the three-qubits 
XXZ-model |g and the XY-model \l% . Further top- 
ics are ent ang lement close to quantum phase transitions 
H EI El EH El E3 and global entanglement with an 
application to quantum error correction code subspaces 

In the present paper, several new aspects of quantum 
entanglement are discussed, in particular how the var- 
ious measures of entanglement can be related to cor- 
relation functions. After introducing briefly the basic 
notations and definitions in the next Section, the func- 
tional dependences of the entanglement measures concur- 
rence 1 1 OL |2CI . i-concurrence [2l| (in small systems) and 
3-tangle | 22| on spin correlation functions (including spin 
expectation values) are established in Sec. lIIII Necessary 
and sufficient conditions for a positive concurrence are 
found. In Sec. IIV1 the expectation values, correlation 
functions and concurrence of both ground and excited 
states of the one-dimensional XXZ-model as well as the 
mixed state of the quantum system at finite temperature 
are calculated analytically in terms of the eigenenergies. 
The concurrence of a N = 4 quantum spin system and 
the critical temperature where the concurrence vanishes 



are examined in detail. Results are also presented for 
N = 2, 3, 5 and 6 qubit systems. Finally, the entan- 
glement of a quantum system with N — 5 qubits in a 
NMR quantum error correction experiment ^23j is dis- 
cussed and partly quantified in terms of the entanglement 
measures in Sec. W\ 



II. BASIC NOTATIONS 

Consider a quantum system consisting of N qubits on 
numbered sites. The basis of the state of one qubit is 
given by |0), |1) which are the eigenstates of a z (<j x , 
a y , a z denote the Pauli spin operators) with eigenval- 
ues — 1, +1, respectively. An unentangled state of N 
qubits is the direct product of the single qubits, e.g., 
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unambiguous then indices indicating site numbers will be 
omitted in the following because the qubits are arranged 
with increasing site number. Thus site information is 
contained in the ordering of the qubits. The Hamilto- 
nian H and the density operator p describing such quan- 
tum spin systems are usually expressed in terms of the 
identity operator /, the Pauli spin operators, and/or the 
operators :— | (a x ± ia v ). 

The state of the spin system becomes mixed at fi- 
nite temperatures. The operator representing this state 
is frequently called thermal density operator. In ther- 
modynamical equilibrium, it is given by the operator 
p = Z^ 1 exp[— /3-H], where f3 = (AisT) -1 , Ub denotes the 
Boltzmann constant, T is the temperature of the system 
and Z — Trexp[— f3H] is the partition function. 

Spin expectation values and correlation functions are 
defined as 
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where n, 



=Tr(p<---<) 
G {!,..., N} and v, . . 



(1) 



2 



{x,y, z,+,— } specify qubit and operator, respectively. 
Furthermore, in what follows, the z-component of the 
total spin operator S z :— J2n=i <J rn the spinflip oper- 
ator F := <3) n= i G n anc i assuming periodic boundary 
conditions, the translation operator T(l) := exp[— ilaK] 
with the lattice constant a and the momentum operator 
hK = P will be used occasionally. 



For a density operator with the coefficients ©, one 
has 



Ai = A 2 = ^ 
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III. ENTANGLEMENT AND CORRELATION 
FUNCTIONS 

The functional dependence of entanglement (measured 
in terms of the concurrence, i-concurrence and 3-tangle) 
on correlation functions of the operators a x , a v , a z , ct ± 
is now discussed as far as possible without an explicit 
specification of the model Hamiltonian. 

Using the basis |0) and |1), the expansion coefficients 
of the (reduced) density operator of one qubit n (1 < n < 
N) are given by spin expectation values only: 
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In the same manner, the (reduced) density operator of 
two qubits n and 77i(l<rt<TO<iV) can be expressed 
in the basis |00), |01), |10) and |11). If the Hamiltonian 
commutes with the z-component of the total spin op- 
erator, the corresponding expressions can be simplified, 
yielding 
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and all other coefficients are equal zero. 

Concurrence C has been introduced by Wootters pol | 
as a measure to quantify entanglement. Let p be the 
density operator representing a pure or mixed state of 
two qubits n and m. Then 
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Cnm — 2A maa: ^ ^ Aj , 
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where X max := max (Ai, A 2 , A3, A 4 ) and Ai, A 2 , A3, A 4 
are the non-negative, real eigenvalues of the matrix R = 
y/p (o-v ® o-y) p* (ay ® &u). 
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Thus Eqs. and JJJ) yield the functional dependence 
of the concurrence on correlation functions using S z - 
symmetry only. 

Cases 3 and 4 of Eq. (J7J are not interesting because 
C n m < and thus C nm = 0. With the help of cases 1 
and 2, it is straightforward to find the following necessary 
and sufficient conditions for entanglement, 
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respectively. These results are similar to the conjecture 
that the ground state of the transverse Ising-model and 
the XY-model is entangled, iff, according to [14|, Kf^ — 

Eqs. © and © can be interpreted in the following 
way: If the state of two qubits in a system with K% — 
and/or = is entangled then the z-components 
of the spins must be correlated antiferromagnetically. 
The maximal entangled states are the two Bell-states 
|^±) = ^ (|01) ± 1 10)). If K z n K z m > 0, e.g., if an appro- 
priate external magnetic field is applied, entanglement 
of qubits with ferromagnetically correlated z-components 
of the spins is possible. The sufficient condition requires 
moreover that the correlations of the two qubits need to 
be greater than a minimum value to create entanglement. 
Again an appropriate external magnetic field reduces this 
demand. 

If the system exhibits additional spin flip symmetry, 
K z = K z m = and K+ m = K~+ result. Then Eqs. ©, 
© and J7J simplify and case 1 of Eq. (7J coincides to 
the result published in Necessary and sufficient con- 
ditions for entanglement are now 
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respectively. Here the relation Kf^ = K^ m — 2K+ n , 
which is correct because of S z and F symmetry, was used. 

I- concurrence C has been proposed by Rungta et 
al. |2l| as an entanglement measure. Let AB be a quan- 
tum system consisting of two subsystems A and B with 
dimensions oIa and ds, respectively. The density opera- 
tors representing the state of these systems are denoted 
Pab, PA and pb, respectively. If pab represents a pure 
state then the entanglement of this state with respect to 
the two subsystems A and B is quantified by 



IV. XXZ-MODEL 
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where pa — Tr^ (pab) is the reduced density operator 
of subsystem A. It is known from [2l| that < Ca-b < 

\J 2^f-, where d = mux (dA, ds)- A different notation is 

occasionally used for qubits: For example C12-34 denotes 
the entanglement of the state where subsystems A and 
B consist of qubits 1, 2 and 3, 4, respectively. Note that 
C nm — C n ~m if the state of qubits n and m is pure. 
From Eqs. J5J) and (fT2*)l . it follows that 



Cn-rest = \f 1 - (K Z n ) 2 - AK+ Kn 
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If the Hamiltonian commutes with S z , Eqs. and 112(1 
yield 
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In an analogous way the i-concurrence of three and more 
qubits can be expressed in terms of correlation functions. 

Two highly entangled qubits cannot be much entan- 
gled with the remaining system and vice versa. This 
property is ensured in Eqs. (fP3)) and ((Fiji . They indicate 
high entanglement in the system if the absolute values of 
expectation values and correlation functions are as small 
as possible (preferable zero). This is contrary to the re- 
quirements for a high concurrence. 

3-tangle r has been suggested by Coffman et al. [22| to 
quantify the entanglement of a pure state of three qubits 
1, 2 and 3 in the following way: 



T 123 — Cl-23 — C\ 2 — C\ 



(15) 



where Cf_ 23 = 4det(pi) = Cf_ 23 and p x = Tr 23 (pi23)- 
Note that T123 does not contain the entanglement of two 
out of the three qubits and T123 does not depend on the 
arbitrary choice of qubit 1 as the "central" qubit. 

The 3-tangle T123 can be expressed in terms of correla- 
tion functions if the Hamiltonian of the system commutes 
with S z . This is achieved by expressing the right hand 
side of Eq. (|15|) in terms of correlation functions with the 
help of Eqs. Q, and (H2J). 



The Hamiltonian H(J, A) of the one-dimensional (spa- 
tial) homogeneous XXZ-model reads (cf. (2^|) 
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The coupling constant J specifies the strength of nearest- 
neighbor spin interaction. Anisotropy in spin space is 
quantified by A. Periodic boundary conditions are as- 
sumed. In what follows, all energies are measured in 
units of J. 

The XXZ-model possesses some interesting symme- 
tries. The Hamiltonian (|16fl commutes with the z- 
component of the total spin operator S z , the spinflip 
operator F and the translation operator T(l). Unfortu- 
nately S z and F do not commute but of course it is possi- 
ble to classify eigenstates of H by eigenvalues s of S z and 
eigenvalues k of ^ln[T(l)]. Because of F-symmetry, it 
is sufficient to solve the eigenvalue problem of H in sub- 
space with s < 0. 

It was shown in [2{| that H{ J, A) and H(—J, —A) pos- 
sess for even N a spectrum of identical eigenvalues in each 
subspace of s because the operator A := &)^=13 a n 
commutes with S z and AH(J,A)A^ 1 = H(-J, - A) = 

//•:•/• A;. 

Some correlation functions of the XXZ model are in- 
terdependent. If only eigenstates with equal s participate 
in the thermal density operator then it is straightforward 
to show that 
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, m 2 are the elements of 



N, Mi U M 2 



Mi and M 2 , respectively, £1 + £2 
{1,2,..., N} and M x n M 2 = 0. 

If H has S z - and F-symmetry, only K^ m and 
appear in Eq. (JJJ . These correlation functions can be ex- 
pressed in terms of the partition function. For example, 
Kl+i) ^d K+~ +1) read (cf. 0) 
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Using these relations, the correlation functions and con- 
currence of the eigenstates and the thermal state of near- 
est neighbor qubits can be calculated by knowing only the 
eigenvalues of the Hamiltonian. It is straightforward to 
express further expectation values and correlation func- 
tions in terms of the partition function using the same 
method. Possibly, the Hamiltonian has to be supple- 
mented (e.g. adding to H appropriate external magnetic 
field terms yields K z again as derivatives of InZ). 

As another application of eq. 0, the concurrence 
of nearest neighbor qubits of the ground state in the 
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anisotropic XXZ-model with J = — I, A — — | and an 
odd number of qubits is considered. It is known from p7^ 

— 2 "I- 2N 2 dna ^„ (rl+1) - 2 1X n(n+l) ~ 

T6 + T6lW- Therefore ATf ?„_,_,•, < and 



that A^„ +1) 

2 n(n+l) " 



1 < 

the concurrence is C 
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L n(ra+1) 

for TV > 3. Thus 
) . Concurrence is 



_ 3 (-i _ 1 
n(n+i) — g W 1 
increasing with odd N whereas the concurrence of near- 
est neighbor qubits of the ground state in the isotropic 
antifcrromagnetic Heisenberg model decreases with in- 
creasing even N in all cases that have been calculated by 
O'Connor et al. j|. 

Now the XXZ model is considered on a finite chain. 



Of course, the calculation of eigenstates and eigenvalues 
is getting more involved with increasing N in general. 
Therefore, in what follows, only small spin chains with 
2 < N < 6 are considered. 

For the case N = 4, the eigenstates \4>) are given in Ta- 
ble H] together with C„(„ +1 ) and C„(„ +2 ), i.e., the entan- 
glement of nearest and next-to-nearest neighbor qubits in 
these eigenstates measured in terms of concurrence (@J. 
Eigenstates with s > are obtained by applying F on 
eigenstates with s < 0. 

The partition function, correlation functions and con- 
currences at finite temperatures are calculated as 
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where C := e' 3 ' 7 and /z 1)2 := -±A =p ^VA 2 + 8. 

The concurrence C„( n+ i) of the state of two nearest 
neighbor qubits as a function of anisotropy A and tem- 
perature T is depicted in Fig. ^ The energies together 
with the concurrences of the individual eigenstates are 
responsible for all described features. At T = 0, the 
change of the ground state from E = A (s = ±2, 
k = 0) to E = \i\ (s = 0, k = 0) causes the discon- 
tinuity at A = — 1. The position of the maximum in 
C n (n+i){A,T = 0) is at A = A max — 1. With increas- 
ing temperature, A max increases but C n ( n+1 )(A max ,T) 
decreases monotonously. Concurrence C n i n +x) f° r fixed 
A is a monotonously decreasing function of temperature. 
As more energies near the ground state energy exist as 
quicker decreases concurrence with temperature. For ex- 
ample, the plateau region in the dependence of C n („ +1 ) 
on T for A > 4 stems from the with A increasing gap be- 
tween ground state energy and most energies of excited 



states. The critical temperature T c is defined as the low- 
est temperature above which the entanglement measure 
(here the concurrence) indicates an unentangled (part of 
the) state (cf. [5^, p. 155]). It is easily identified as the in- 
tersection of the zero-surface and the surface of the func- 
tion C n („ + i) in Fig.^ The projection of the critical tem- 
perature T c and the lines of equal C n t n +x) are depicted 
in the lower part of Fig. ^ In this way it is easy to iden- 
tify parameter regions of states with a certain minimal 
entanglement. Note that lines of finite equal concurrence 
are not increasing monotonously with increasing A but 
T c does. 

In Fig. |21 the critical temperature T c of the entangle- 
ment (measured in terms of concurrence) of the state of 
two qubits in the XXZ-model ( J ^ 0) for 2 < N < 6 as 
a function of anisotropy A is shown. 

The transformation J — » — J and A — > — A leaves the 
critical temperature invariant for even N. If \ip) is an 
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TABLE I: Classification of the eigenstates of the XXZ-model (N — 4) and concurrence of nearest and next-to-nearest neighbor 
qubits. Normalization factors are given as 771,2 := y4 + 2(^i,2) 2 , where mi,2 := — ^A =p ^r\/A 2 + 8. 
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eigenstate of H (J, A) with eigenvalue E then |</>) = A |t/>) 
is the corresponding eigenstate of H(—J, —A) with the 
same eigenvalue and identical entanglement because A is 
a local unitary transformation and entanglement is in- 
variant under local unitary transformations. Thus the 
thermal density operators of both Hamiltonians are uni- 
tary equivalent and possess identical entanglement and 
critical temperatures. No such symmetry exists for odd 
N. Actually, for the states of nearest neighbor qubits 
(N = 3) and next-to-nearest neighbor qubits (N = 5) 
entanglement is only possible for J < 0. In all consid- 




FIG. 1: The 3D-plot shows the concurrence C n ( n+ i) of the 
state of two nearest neighbor qubits in the XXZ-model (N = 
4, J > 0) as a function of anisotropy A and temperature T. 
The 2D-plot shows the projection of critical temperature T c 
( ) and lines of equal C n ( n+ i) ( ). 



ered cases the inequality T C (N, J < 0) > T C (N, J > 0) is 
valid. 

One observes in Fig.Hthat T c = for (J/ \J\)A < -1 
independently of N and the choice of the two qubits. It is 
known from [24| that for all N, J < and A > 1 the two 
eigenstates of the Hamiltonian (|16|) with s = ±A are 
ground states. These ground states are not entangled 
and they cause the thermal state to be unentangled for 
all temperatures. The same reasoning applies for even 
N, J > and A < —1, because of the symmetries of 
the XXZ-model with periodic boundary conditions. The 
ground state may change at different A for odd N and 
J > (e.g. at A w -0.809015 considering the XXZ- 
model with N = 5 and J > 0). 

Furthermore critical temperature of geometrically 
equivalent aligned qubits is decreasing with increasing N 
for even N. This tendency is consistent with the depen- 
dence of concurrence on N in the isotropic Heisenberg 
model with an applied external magnetic field (cf. Q). 



V. ANALYSIS OF AN EXPERIMENT 

Finally, the entanglement of the state of the quan- 
tum system in a NMR-experiment about quantum er- 
ror correction |23| is quantified in terms of concurrence, 
i-concurrence and 3-tangle. Five qubits are provided by 
different atoms in 13 C labeled transcrotonic acid (synthe- 
sis and properties, see |2|j) solved in deuterated acetone. 

One molecule can be approximately described by the 
one-dimensional spatial inhomogeneous XXZ-model in- 
cluding an external magnetic field because the coupling 
constants of non-neighboring qubits are much smaller 
than the coupling constants of nearest neighbor qubits 
(see [H i^). The Hamiltonian H(J n ,A,u n ) of this 
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(J/\J\) A 



FIG. 2: Critical temperature T c of the concurrence of the 
state of two qubits in the XXZ-model ( J ^ 0) for 2 < N < 6 
as a function of anisotropy A. Panel a) shows nearest neigh- 
bor qubits for N = 2 ( ), N = 4 ( ) and AT = 6 

( ); next-to- nearest neighbor qubits for N = 4 ( ) 

and N = 6 ( ); next-to- next-to- nearest neighbor qubits 

for N = 6 ( ). Panel b) displays nearest neighbor qubits 

for TV = 3 (J < 0: ; J > 0: identical zero) and N = 5 

(J < 0: ; J > 0: ); next-to-nearest neighbor qubits 

for N = 5 (J < 0: ; J > 0: identical zero). The in- 
sets give the dependence of these functions at larger values of 
(J/ | J|)A. Of course, the entanglement vanishes in the Ising 
model limit of 1161 . i.e., for |A| — » oo. 



model reads 

1 5 

11=1 

where the coupling constants J n (n = 1,...,4) specify 
the inhomogeneous strength of nearest neighbor inter- 
action, A determines the anisotropy in spin space and 
the effect of the external magnetic field is included in 
uj n = ujP + oj^ (n = 1, . . . , 5) which are the sums of pre- 
cession frequencies cu? and chemical shifts for each 
individual qubit (data in [ISIi^l)- Of course, now open 
boundary conditions are applied. 

The five-qubit code for quantum error correction is 



used to encode qubit 2 in the experiment. The encoding 
is shown in Fig. |3| The quantum system is in a highly 
mixed state, i.e., the coefficients of the density opera- 
tor are close to the coefficients of the identity operator, 
because the experiment is performed at room tempera- 
ture. In the beginning, the quantum system is prepared 
in a way that only molecules in the initial state 1 11111) 
give a signal on NMR measurements. Then one says that 
the quantum system is in the pseudo-pure state 1 11111) 
(Ref. 30]). The pseudo-pure state 1 11111) is an eigen- 
state of the Hamiltonian H27|) as well as the Hamiltonian 
including all interactions of qubits and the applied exter- 
nal magnetic field described in [2^, [2|j . Furthermore it 
is an cigenstate of S z . Thus going to a frame of refer- 
ence that rotates around the z-axis does not change the 
density operator of the initial state (see 0, p. 287]). 

The pseudo-pure state of the quantum system at sev- 
eral stages (A, B, C, D and E, cf. Fig. |3J| during en- 
coding was calculated by the product-operator-formalism 
(see |3l|, chapter 11]). Therefore, the conservation of 
the pseudo-purity of the state of the quantum system 
is assumed, i.e., there is no interaction between different 
molecules and encoding is implemented so quickly that no 
decoherence occurs. The results are given in Table ITTlto- 
gether with the expectation values K* and if+ = (K~) 
(with n = 1, 2, 5). 

It is straightforward to calculate the entanglement of 
one qubit with the remaining qubits by inserting these 
expectation values into Eq. (|13f) . In this way it is easy to 
get a quick overview about the possible entanglement in 
the quantum system. Note that it is not appropriate to 
use Eqs. Q and Q or (|14"}> here because the pseudo-pure 
state does not comply with the necessary S^-symmetry 
in general. 

The pseudo-pure state at the various stages is now dis- 
cussed in detail: The initial state is not entangled. At 
position A, the state is not entangled as well. So far only 
local operations have been performed and these cannot 
create entanglement. 

At position B, qubits 1, 2 and 5 are not entangled but 
C34 = 1 . Actually, the state of qubits 3 and 4 at position 
Breads |^) 34 = \ (|H) + 1 10> - |01) + |00>) and it is con- 

encoding 



A B CD E 




FIG. 3: Encoding of qubit 2 based on the five-qubit code. The 
horizontal lines represent the qubits. The gates denoted 8a" 
and 180° cr^of implement e _ 5 9<T ? and e~i na ^ a ^ , respectively. 
Here a £ {x, y, z} and a, b £ {1,2, 3, 4, 5}. 
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TABLE II: Pseudo-pure state \ip) of the quantum system at several states during encoding. The expectation values K„ 
and Kt = (Kn)* (with n = 1, 2, . . ., 5) are given for each state. Notation: |l a > := 4- (|1) + |0)), \0 X ) := 4j (|1) - |0>), 

:= rfc (|1> + » |0», |0y) := 4j (|1> - * |0», |1.) := |1) and |0.) := |0). 



Position 


w 






Kg 


Kt 


Kl 


Kt 


Kt 


Kt 


Kt 


Kt 


A 







1 








1 


i 

2 





2 


~ 2 





B 


^|lyl*)®(|l a l») - 10,0,=)) 





1 








1 


2 














C 


^75 |ly) ® (10*1*0*) + llzOzlz}) ® |l z ) 














1 


2 














D 


| |l y > ® [\Q*1 X ) ® (|l»0x) + 10,1*)) - i|l*0„) ® (|1»0 S ) - |0*1*))] 

















2 














E 


5^5 fll*) ® [|0 V 1») <8> (|1*1*> + j0 z z )) - \l y z ) (8> (|1*1*) - |0«0 Z »] + 
i 10,) <g> [|0j,0 e ) <8> (|1 Z Z ) + jO z l z )) - \l v l x ) ® (|1 Z Z ) - |0 Z 1 Z ))]} 

































form to the Bell-states 1^) and 1^) = -4= (|00) ± 1 11) ) 
up to a local unitary transformation. 

At position C, only qubits 2, 3 and 4 are en- 
tangled: C 2 -34 = C3-24 = C4-23 = T 23 4 = 
1, where the state of these qubits reads |t/j) 234 = 
\ (|110) + 1 101) - |010) + |001>). It is conform to the 
cat-state ^= ( 1 000) + |111)) up to a local unitary trans- 
formation. Two out of these three qubits are not entan- 
gled as usual for a cat-state. 

At position D, only qubit 1 is not entangled. 
The state of the remaining qubits is conform to 
\ (|0110) + |0101) - i |1010) + i |1001>) up to a local uni- 
tary transformation. The analysis of qubits 2, 3, 4 and 5 
shows no entanglement of the state of two of these qubits. 
The entanglement of a state of three qubits cannot be 
calculated because tracing off a qubit generates in gen- 
eral a mixed state and i-concurrence can only be applied 
to pure states. But it is C2-345 = C3-245 = C4-235 = 

C5-234 = 1, C23-45 = 1 and 624-35 = ^25-34 = 

At the end of the encoding sequence (position E), 
all qubits are entangled: Ca-b = 1 if A indicates 
one arbitrary qubit and B the remaining four qubits; 

Ca-b — yf if A indicates two arbitrary qubits and B 
the remaining three qubits. Again there is no entangle- 
ment of the state of two qubits and the entanglement 
of a state of three or four qubits cannot be quantified 
so far. These results coincide with the ones in 0|. It 
was already pointed out there that all states in a certain 
fife-qubit error correction code subspace possess maximal 
global entanglement but vanishing concurrences. 

Clearly, in this experiment, entanglement is created 
during encoding and it expands in a geometrical sense, 
i.e., the number of qubits involved in the entanglement 
increases with the progressing encoding sequence. 

Unfortunately, it is not possible to quantify the entan- 
glement of the state at positions D and E completely be- 
cause of the lack of suitable measures. But all calculated 
i-concurrences exhibit their maximal values at position 



E. Thus it is a reasonable conjecture that an entangle- 
ment of four or less qubits does not exist there because 
entanglement cannot be shared arbitrarily (cf. p^ h 



VI. SUMMARY 

The entanglement measures concurrence, i- 
concurrence (for one or two qubits in one subsystem) 
and 3-tangle have been successfully expressed in terms 
of correlation functions. In addition, necessary and 
sufficient conditions for a positive concurrence have 
been formulated. These results have been used in the 
remaining paper because they can simplify calculations: 
The concurrence of eigenstates or the thermal state 
have been calculated analytically knowing only the 
energies of the eigenstates and their dependences on 
the parameters of the system. Furthermore potential 
quantum entanglement in a quantum system has been 
detected by the examination of spin expectation values. 

A detailed analysis of concurrence and critical temper- 
ature in the XXZ-model with 2 < N < 6 qubits has been 
accomplished. 

Finally, the entanglement of the state in a NMR- 
experiment has been discussed quantitatively. Different 
kinds of entanglement have been identified. This calcu- 
lation shows the relevance of entanglement measures in 
actual experiments because they allow an analysis of the 
importance of entanglement for the quantum-algorithms. 
Despite the information, which is obtained with the avail- 
able measures, further measures are needed for a com- 
plete insight. 

The entanglement measures might be useful design- 
ing new experiments (possibly utilizing advanced types 
of qubits, e.g., spin cluster qubits [32^ ) that set up states 
with different entanglement and prove or disprove the 
benefit of entanglement in different quantum-algorithms. 

One of us (H.F.) thanks John Schliemann for useful 
discussions. 
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